pub fn new_birthday_probability(n: u32) -> f64 {
    let mut mul = 1f64;
    for i in 0..n {
        mul *= (365f64 - (i as f64)) / 365f64;
    }
    1f64 - mul
}
